QCD Factorization for Deep-Inelastic Scattering 
At Large Bjorken x B ~ 1 - (9(A QC d/Q) 



Panying Chen, 1 Ahmad Idilbi, 1 and Xiangdong Ji 1,2 

1 Department of Physics, University of Maryland, 
College Park, Maryland 20742, USA 
2 Department of Physics, Peking University, Beijing, 100871, P. R. China 

(Dated: February 2, 2008) 

Abstract 

We study deep-inelastic scattering factorization on a nucleon in the end-point regime xb ~ 
1 — C(Aqcd/Q) where the traditional operator product expansion is supposed to fail. We argue, 
nevertheless, that the standard result holds to leading order in 1 — xb due to the absence of the scale 
dependence on (1 — xb)Q- Refactorization of the scale (1 — xb)Q 2 in the coefficient function can be 
made in the soft-collinear effective theory and remains valid in the end-point regime. On the other 
hand, the traditional refactorization approach introduces the spurious scale (1 — xb)Q in various 
factors, which drives them nonperturbative in the region of our interest. We show how to improve 
the situation by introducing a rapidity cut-off scheme, and how to recover the effective theory 
refactorization by choosing appropriately the cut-off parameter. Through a one-loop calculation, 
we demonstrate explicitly that the proper soft subtractions must be made in the collinear matrix 
elements to avoid double counting. 



I. INTRODUCTION 



In lepton-nucleon deep- inelastic scattering (DIS), the Bjorken regime with virtual photon 
mass Q 2 — > oo and xb = Q 2 /2Mu fixed presents a textbook example of perturbative QCD 
(pQCD) factorization p|. In this regime, the scale (1 — x) a Q goes to infinity (we drop the 
subscript B on x henceforth), where a is any real number, or at least much larger than 
the soft QCD scale Aq CD . An alternative DIS regime is Q 2 — ► oo with (1 — x)Q ~ Aqc D , 
where the final hadron state invariant mass (1 — x)Q 2 ~ QAqcd is still large and is distinct 
from the resonance region. This large- X regime has so far received little attention in theory, 
possibly because it covers only a small kinematic interval in real experiments. The existing 
QCD studies in the literature are somewhat controversial 0,0]. 

In this paper, we present a factorization study of DIS at this regime. The main point we 
advocate here is that the standard pQCD factorization remains valid in this new kinematic 
domain to leading order in 1 — x. Then we move on to discuss refactorization which factorizes 
the physics at scale (1 — x)Q 2 from that at Q 2 . The useful theoretical framework for this 
purpose is the so-called soft-collinear effective field theory (SCET) developed recently 0. 
Indeed, the first treatment of the large-x region in DIS using SCET was made in p] and 
followed in 0,0] for (1— x)Q ~ Aqcd [see also 0.] Because the scale (1 — x)Q does not enter 
in the perturbative calculation, the final result amounts to a standard pQCD factorization, 
with the additional benefit that the refactorization becomes manifest. One subtlety we 
discuss extensively in this paper is the role of the soft contribution and its relation to 
the light-cone parton distribution. In a recent paper 0, a factorization formula for the 
DIS structure function is derived in SCET, which is similar to what we find here. However, 
because of the lack of a consistent regularization and clear separation of contributions among 
different factors, the result does not recover that of Ref. p in the limit (1 — x)Q ^> Aq C d- 

The traditional approach of refactorization was pioneered in Q, where a new parton 
distribution together with soft factor and jet function is introduced. These matrix elements 
are designed to absorb large logarithms generated from soft gluon radiations off on-shell 
lines. Evolution equations for them are derived and solved to resum the large soft-gluon 
logarithms. In this approach, the factorization of scales are not apparent from the start. 
Moreover, a spurious scale (1 — x)Q appears in all factors which makes them nonperturbative 
in the regime of our interest. After reviewing this, we present a more general factorization 
along this line with dependence on a rapidity cut-off p. Different choices of the cutoff lead 
to redistributions of large logarithms in different matrix elements. A particular choice yields 
a picture similar to that of the SCET approach. 

The presentation of this paper is as follows. In section II, we argue that the standard 
factorization approach remains valid in the end-point regime x ~ 1 — 0(Aqcb/Q)- In section 
III, we present an effective field theory approach to refactorization, following the previous 
work of Ref. 0- The difference is that our result is consistent with that of Ref. 5|, with 
the jet factor absorbing all physics at the intermediate scale (1 — x)Q 2 . We explain that 
the soft and collinear contributions combine to give the light-cone parton distribution. In 
section IV, we first review the traditional factorization in which various matrix elements 
are introduced to account for soft gluon radiations. We then show how to derive a more 
general factorization with a rapidity cutoff. We demonstrate explicitly that the proper soft 
subtractions must be made in the collinear matrix elements to avoid double counting. By 
choosing different cut-off, we find different pictures of factorization and large-logarithmic 
resummation. The effective field theory refactorization can be recovered this way. 
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II. VALIDITY OF THE STANDARD QCD FACTORIZATION AT xb ~ 1 — 

0(A QCD /Q) 



The standard pQCD factorization theorem is derived in the Bjorken limit in which Q 2 — > 
oo and a; is a fixed constant between and 1. To leading order in l/Q 2 , the proton's 
spin-independent structure function Fi(x,Q 2 ) can be factorized as 

F 1 (x,Q 2 ) = j:[ 1 ^C f (^9p jqf (y,, 2 ), (1) 

where /i is a factorization scale, Cf is the coefficient function depending on scale Q 2 and 
fi 2 (factorization scale), and qf is a quark distribution of flavor /. For simplicity, we omit 
the quark charges and gluon contribution which are inessential for our discussion. In the 
moment space, the factorization takes a simple product form, 

f 1 n (Q 2 ) = J2 c MQ 2 /^M^ > ( 2 ) 

/ 

where the moments are defined the usual way, e.g., q^ = dxx N ~ 1 (q(x) + q(x)) with N 
even. 

The above factorization in principle is invalid in non- Bjorken regions. However, we argue 
that it still holds in the regime of our interest, (1 — x)Q ~ Aqcd, with the same physical 
parton distributions to leading order in 1 — x. The main point is that although we now have 
a new infrared scale (1 — x)Q, it does not appear in the above factorization. 

Indeed, as x — * 1 the coefficient function Cf can only depend on the two hard scales — 
invariant photon mass Q 2 and the final-hadron-state invariant mass (1 — x)Q 2 , both remain 
large when (1 — x)Q ~ Aq CD . Hence there is no emerging infrared scale entering physical 
observables in this new regime, and the original factorization remains valid to leading order 
in 1 — x. 

The above observation can be seen clearly in the one-loop result. The coefficient function 
at order a s in the x-space 



2\n((l-x)Q 2 /i2 2 ) -3/2\ /3 Q 2 9 



l-x J, \2 fj, 2 2 3 
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ln-% + - + — U(l -x) 



, (3) 



where we have neglected higher order in (1 — x). The scheme we use here is the modified 
minimal subtraction (MS). In term of moments, one finds 



r (i)/Q 2 Q 2 \_a 



Nfi 2 2 Nfi 2 fi 2 fi 2 2 6 



(4) 



The scale dependence is manifest: The first two logarithms come from physics at scale 
[i 2 = Q 2 /N, whereas the next two logarithms come from scale Q 2 . Clearly, there is no 
physics from scale /i 2 ~ (1 — x) 2 Q 2 ~ Q 2 /N 2 . Therefore, even when (1 — x)Q becomes of 
order Aqcd, the coefficient function has no infrared sensitivity to it. 

The fundamental reason for the absence of the scale (1 — x)Q in a physical observable is 
that it is not a Lorentz scalar, whereas (1 — x)Q 2 is the invariant mass of the final hadron 
state. In principle, the energy of soft gluons and quarks in the Breit frame is of order 
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(1 — x)Q which can appear in the factorization. However, this happens only for frame- 
dependent factorization. The factorization we quoted above is frame-independent and thus 
any non-Lorentz scalar cannot appear. 

Although the above conclusion appears simple and natural, we have not seen it stated 
explicitly in the literature. 



III. REFACTORIZATION: EFFECTIVE THEORY APPROACH 

In the large-x region, independent of whether (1 — x)Q ^> Aqcd or ~ Aqcd, there is an 
emerging "infrared" scale (1 — x)Q 2 <C Q 2 . Of course, we always assume (1 — x)Q 2 ^> Aq CD . 
The presence of this new scale suggests a further factorization in which the physics associated 
with scales Q 2 and (1 — x)Q 2 is disentangled. This type of factorization was proposed by 
G. Sterman and others for the purpose of summing over the large double logarithms of 
type aj[ln*(l — x)/(l — x)] + (i < 2k — 1) in the coefficient functions Q. We consider this 
refactorization in this and the following sections, commenting on its applicability in the 
region of our interest. 

We first study refactorization in the effective field theory (EFT) approach in this sec- 
tion, and will discuss a more intuitive approach in the next. The EFT method is based on 
strict scale separation, very much like the usual QCD factorization discussed in the previous 
section. When the scales are separated, one can sum over large logarithms by using renor- 
malization group evolutions between scales. Some of the basic discussions here follow Refs. 



EFT Refactorization 



To understand the EFT factorization heuristically, we write the one-loop coefficient func- 
tion in Eq. (@J) in a factorized form, 



r a)(Q 2 Q 2 
N V'aV 



N \ a r 2 



(5) 



where is A-independent and comes from physics at scale Q 2 and N = Ne lE where 7^ 
is Euler constant. The one- loop result for C is 



Q 1 
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(6) 



where the constant term is, in principle, arbitrary; we choose it to be consistent with the 
effective current below. The two-loop result for C can be found in [ljj • The other factor 
comes entirely from physics at scale (1 — x)Q 2 ~ Q 2 /N, 
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(7) 



and the second order result for A^tv can also be found in The key point is that the 
above refactorization of scales works to all orders in perturbation theory and EFT provides 
a formal approach to establish this: The physics at scale Q 2 can be included entirely in \C\ 2 
and that at the other scale is in M. N . 
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To arrive at the above refactorization, we start off at the scale Q 2 at which perturbative 
physics involves virtual gluon corrections to the hard interaction photon vertex. Note that 
the soft-gluon radiations off a hard vertex are usually high-order effects in 1/Q 2 and can 
be neglected. Thus the physics at Q 2 can be found from just the quark electromagnetic 
form factor. Integrating out physics at scale Q 2 is equivalent to matching the full QCD 
electromagnetic current to an effective one involving just soft-collinear physics. 

Jfa = C(Q 2 /^ S (Q 2 /^ 2 ) , (8) 

with the one-loop result given in Eq. JBJ). 

We can run the effective current from scale Q 2 to scale (1— x)Q 2 using the renormalization 
group equation 

/ Jeff( ^ 2/V) = -7i(a.(/0) J cS (Q 2 /^ 2 ) , (9) 

where the anomalous dimension can be calculated from C, 71 = fidlnC/dfi, and has the 
following generic form, 

71 = A(a s )\nQ 2 /fj, 2 + B 1 (a s ) , (10) 

in which A and B\ are a series in strong coupling constant a s and are now known up to 
three loops . 

At scale \x\ = (1 — x)Q 2 , we follow Ref. matching products of the effective currents to 
a product of the jet function, collinear parton contribution, and soft distribution in SCET. 
Introducing a small expansion parameter A, with X 2 Q ~ Aqcd, 1 — x in the region of our 
interest scales like A 2 . The collinear partons at the matching scale /1 2 - have momentum 
(p + ,p~p±) ~ Q(l, A 4 , A 2 ) [our notation for light-cone components for arbitrary four-vector / 
is I = (/ + , l~ , l±) with = -j^ip ± ^ 3 )], and the soft partons have momentum (A 2 , A 2 , \ 2 )Q. 
The moment of the structure function F\ after the second stage matching has the following 
form 0|, 

if (Q) = C 2 (Q 2 /ri)J P (N, Q 2 //j? I )(j) N (ii 2 I )S N (iJ l 2 I ) , (11) 

where various factors are defined as follows. 

The jet function Jp(N, /i 2 ) is related to the absorbtive part of the hard collinear quark 
propagator Q P , 

(0|T [W^ n (z)UW n (0)} |0) = i-*=J ^e- ikz g P (k) . (12) 

where ^ n is a collinear quark field and W n is a Wilson line along the light-cone direction n 
(n 2 = n 2 = 0, n-n = 1). A hard collinear quark has momentum p+k, where p is the so-called 
label momentum with p + ~ Q and k is a hard residual momentum with components of order 
Q(l, A 2 , A). Therefore, the virtuality of the hard-collinear quark is 2p ■ k ~ \ 2 Q 2 , consistent 
with that of the hadron final state. The jet function has no infrared divergences because the 
hadron final states are summed over. However, it does have light-cone divergences which 
are handled by the standard minimal subtraction method. An important feature of the jet 
function is that it is only sensitive to physics at scale /if. In fact, at one-loop order, the jet 
function reproduces the Ai n function in the previous section. 

The soft contribution in Eq. (12) is defined in terms of the soft Wilson lines 2| : Y n (x) = 



Vexp[—ig f* dsn ■ A us ] and Y n (x) = V exp[—ig dsn ■ A us ], where A us are the so-called 
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ultra-soft gluons with momentum Q(X 2 , A 2 , A 2 ) and V stands for path ordering [the sign 
convention for the gauge coupling is = + ig A^} , 



Tr 



y- T y s 

1 n 1 n u 



X + 



n 



id 



n ■ p 



1 n 1 7 







(13) 



where the ratio in the delta function fixes the momentum of the emitted gluon to be soft 
(of order A 2 ). As such, the soft factor is a non-perturbative contribution. 
The collinear contribution, 



P 



x 



n-V + 
n ■ p 



P 



(14) 



comes from collinear quarks and gluons with momentum (1, A 2 , X)Q and n ■ V + is the total 
light-cone momentum carried by the partons. In Ref. [2j], the collinear contribution was 
identified as the usual Feynman parton distribution. This, however, is incorrect because the 
soft gluons with longitudinal momentum (1 — x)Q ~ A 2 in the proton cannot be included 
in the collinear contribution according to the definition of the collinear gluons in SCET. On 
the other hand, the Feynman parton distribution contains a factorizable soft contribution in 
the limit x — > 1 k! 3 • Further discussion on this issue will be made in the next subsection 
as well as in the next section. The correct approach is to combine the soft and collinear 
contributions together to get the correct Feynman parton distribution in the x — > 1 limit. 

Thus EFT arguments lead finally to the following refactorization, valid when (1 — x)Q ~ 
Aqcd at leading order in 1 — x, 

Fi{Q 2 ) = C 2 {Q 2 /ti)J P {Q 2 /Nti)q N {ti) + Oil - x), (15) 

where gAr(/x 2 ) is the moment of the quark distribution, and the jet function is exactly the 
M.N function introduced in the previous section. Although formally this factorization is 
made at n 2 , the product of factors is independent of it. The claim of non-factorizability of 
DIS in this very regime in Ref. [3] was criticized in . On the other hand, the above result 
seems consistent with that of Ref. if used in the same regime. 

The above factorization allows us to resum over large logarithms. Since the physical 
structure function is \i- independent, we can take /ij in the above expression to whatever 
value we choose. For example, if one sets /xf = Q 2 , all large logarithms are now included in 
the jet function. One can derive a renormalization group equation for Jp Solving this 
equation leads to a resummation of large logarithms. 

Alternatively, with the original scale fi 2 , there are large logarithms in C, which can be 
resummed using the renormalization group equation and the anomalous dimension 71. The 
resulting exponential evolution can be regarded as the evolution of the jet function from scale 
Q 2 to fij. The parton distribution gjv (/•*/) runs from \jt\ to a certain factorization scale \xp 



using the DGLAP evolution [15|]. This running generates the logarithms from initial-state 
parton radiations. 

In the above refactorization, no scale (1 — x)Q appears explicitly although soft and 
collinear gluons in SCET do have reference to that scale. This explains that the factorization 
holds in the region of our interest, namely, when (1 — x)Q ~ Aqcd- 



B. Collinear Contribution in SCET and Double Counting 

SCET is an operator approach designed to take into account contributions from different 
regions in Feynman integrals. Calculations in SCET are sometimes formal if without a 
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careful definition of regulators for individual contributions. Occasionally, the regulators in 
different parts must be defined consistently to obtain the correct answer. Otherwise, one 
can easily lead to double counting. The same issue has been discussed recently in Ref . [TtJ • 
To see the need of consistent regulators in SCET, let us consider the usual quark distri- 
bution in the proton. In the full QCD, the quark distribution is defined as 

q(x) = - [ —e iXx {P\lp(Xn) ^e~ i9 ^ dX ' n - A ( x ' n ^(0)\P) , (16) 
2 J 2n 

where ip{x) and A(Xn) are full QCD quark and gluon fields (here we use the vector n 
with mass-dimension 1). Now suppose the nucleon is moving with a high momentum Q 
in the z direction. The quarks and gluons in the proton, in general, have large k + , and 
small momentum k~ and k±, in the sense that they are collinear to the proton momentum. 
Therefore, one may match the full QCD fields in the above expression to the corresponding 
collinear fields in SCET. 

However, the procedure is incomplete for wee gluons with (1 — x) ~ A/Q. Such a gluon 
has a soft longitudinal momentum and is definitely included in the above gauge link. The 
QCD factorization theorem shows that soft gluons do not make a singular contribution to 
the parton density, but it does not exclude the non-singular wee gluon contributions of type 
[ln fc (l — x)/(l — + . In fact, the wee gluon effect in the x — > 1 limit can be factorized out 
into a soft factor S(x), which is responsible for the large-x behavior of the parton distribution 
[l3L fT3| . Therefore, in SCET it is natural to express q(x) in terms of the product of the soft 
factor and true collinear gluon contribution. 

In Ref. 0, the evolution equation was derived for the collinear contribution 4>(x) and is 
found to be the same as the DGLAP evolution, even in the x — > 1 limit. This could be 
the main motivation to identify 4>(x) as q(x). However, the collinear gluons in the one- loop 
Feynman diagrams can no longer be considered as "collinear" if (1 — x)Q ~ Aqcd, and must 
be subtracted explicitly. This subtraction was not made through certain regulators and 
hence there is a double counting. In fact, once the soft-gluons are subtracted, a collinear 
parton jet shall not have singularity in the limit x — > 1. Likewise, the calculation of the 
soft-factor in Ref. should have a soft transverse-momentum cutoff to include just the 
true soft gluons. Thus, the regulators in the soft and collinear contributions must be made 
consistently to avoid double counting. A consistent scheme of defining the soft and collinear 
contributions for a parton distribution defined with off light-cone gauge link can be found 
in Ref. jl3l ]. We will present another example of a consistent regularization in the following 
section. 



IV. REFACTORIZATION: INTUITIVE APPROACH 

The EFT approach for refactorization is gauge-invariant and all factors are defined at 
separate scales. The resummation is of the simple renormalization-group type. However, 
the physical origin of the large logarithms is not entirely transparent. For example, it is well 
known in QED that the double infrared logarithms are generated from soft radiations from 
jet-like lines. This is not obvious in the EFT approach. 

In the approach introduced by Sterman and others 0], the structure functions are factor- 
ized into different factors which have clear physical significance, although each factor now 
contains multi-scales. Explicit equations can be derived to bridge the scales within these 
factors, which allow one to resum large logarithms. The main shortcoming of this intuitive 
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FIG. 1: The leading reduced diagram contributing to the deep-inelastic structure function in x — > 1 
regime. 

approach is the introduction of a spurious scale (1 — x)Q in each factor, which make them 
nonperturbative in the kinematic region of our interest. 

In this section, we first briefly review Sterman's approach in subsection A. We then 
introduce in subsection B a more general factorization approach along this direction, which 
involves a rapidity cutoff. With an appropriate cutoff, we arrive at a picture similar to that 
of EFT. The example also shows that a consistent soft subtraction must be made to obtain 
a correct factorization. 

A. Sterman's Method 

Consider the lepton-nucleon DIS process in the Breit frame in which the initial and final 
partons have similar momentum but move in opposite directions. In the region x — > 1, 
the final hadron state consists of a high-energy jet plus soft gluon radiations. A so-called 
reduced diagram is shown in Fig. 1, showing the space-time picture of the process. There 
are in principle four different scales which are relevant: virtual photon mass Q 2 , fmal-hadron 
invariant mass (1 — x)Q 2 , the soft parton energy radiated off the proton (l — x)Q, and finally 
the genuine nonperturbative QCD scale A QCD . 

According to the analysis in Ref. 0], the reduced diagram can be factorized into various 
physically intuitive contributions, and the structure function can be expressed as a product 
of a soft factor, a final-state jet function, and a parton distribution, 

f 1 dv P y ~ x dm 
F 1 (x,Q 2 ) = H(Q 2 ) -±<t>(v) -^S(w)J(x,y jW ) . (17) 

Jx y Jo i- — w 

The parton distribution <p{x) is not the usual gauge-invariant one on the light-cone. Rather 
it is defined as 

m = ^J ^e-^(P|g(An) 7 + g(0)|P) , (18) 

in A z = gauge, or equivalently there are gauge links along the z direction going from the 
quark positions to infinity. Because it is not truly gauge- invariant, it is frame- dependent. 
In particular, it can depend on the soft parton energy (1 — x)Q. This parton distribution 
contains contributions of both collinear and soft gluon radiations form the initial state quark, 
thus involving double logarithms. 



S 



Likewise, the jet function is denned as 

d^x 

J ~ Disc 



(2tt)< 



e-»'(0|T[g(*)g(0)]|0> 



(19) 



and normalized to 5(1 — x) at the leading order. Once again the jet is defined in the axial 
gauge and is frame dependent. In particular, it depends on the infrared scale (1 — x)Q as 
well. However, this jet function contains no true infrared divergences. It accounts for the 
collinear and soft radiations from the jet final state. 

Finally, the soft function S is defined as the matrix element of Wilson lines first going 
along the n direction from — oo along the light-cone, then going along n direction to +oo. 
The collinear divergences are regularized by A z = gauge. Therefore, it contains no true 
infrared divergences. 

In the moment space, the refactorization appears in a simple form, 



F((Q) = H(Q 2 /^)4> N (Q 2 /^)S N (Q I /^)J N (Q 1 /^ 



(20) 



and H stands for the hard contribution which comes only from virtual diagrams. Again, 
we emphasize that the factorization follows intuitively from the space-time picture of the 
reduced diagrams. However, one pays a price for this: the breaking of Lorentz invariance 
and introduction of a new scale (1 — x)Q. When this scale becomes of order Aqcd, as this 
is the main interest of the paper, all factors becomes nonperturbative in principle. On the 
other hand, the scale (1 — x)Q is spurious, it should be cancelled out. It is unclear how this 
is achieved at the nonperturbative level. 

The physics of the above factorization is best seen through a one-loop calculation. The 
parton distribution is 
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where P qq (x) is the quark splitting function and 



lnYl - x) 



0,1,2. 



(21) 



(22) 



Apart from the divergent term which is the same as the Feynman parton distribution, there 
are extra constant terms which absorb the soft-collinear gluon contribution. Some of these 
come from the scale (1 — x) 2 Q 2 . In moment space, we have 



4> N (x) 



Q~ 



Q 2 



— (3 - 4 In TV) + In 2 — ^ 2 In — I In 



2 ^ + 21n^--4 



(23) 

One may view the double logarithmic terms as from the initial state radiation. The large 
logarithms resulting from large scale differences can be summed through an x-space evolution 
equation pj. 

The one- loop jet function is explicit finite, 



a. 



J(x) = 5(1 -x) + ^C F 

4:71 
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which involves physics at both scales (1 — x) 2 Q 2 and (1 — x)Q 2 . Both type terms generate 
double logarithms, corresponding to the radiations from the jet. In moment space, the jet 
function becomes, 
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(25) 



Again the large logs can be resummed through an evolution equation for J 0. 
The soft function is also finite and at one-loop; 



S(x) = 5(1 



x) + —C F 



2D, 



\n^)5(l-x) 



(26) 



which contains physics at scale (1 — x) 2 Q 2 , generating a single logarithm. In the moment 
space, it is 

a s „ , Q 2 



Sn = 1 H Cf In 

71 A ll 2 



(27) 



Here the collinear singularity is regulated by gauge fixing. 

When summing over the jet, parton distribution and soft contribution, the soft scale 
(1 — x) 2 Q 2 dependence cancels. We are left with only the physical scale (1 — x)Q 2 . All 
the factors introduced above are sufficient to factor away the singular contributions in the 
structure function at x — > 1 limit. In fact, at one-loop 



F[ l \x,Q 2 ) - (0«(x) + jW(x) + S^(x)) = ^C F 



2tt 



O 2 

-4 + ln^- 
fi 2 j 



5(1 



x 



(28) 



which contains only the 5-function singularity. Therefore, the large double logarithms have 
been absorbed either into the parton distribution or the jet function. This is, in fact, the 
purpose of the intuitive refactorization approach: The double logarithms from the initial 
and final state radiations are made explicit through factorization. 

However, because of the presence of the extra scale (1 — x)Q, the above refactoriza- 
tion is not very useful in the region of our interest because all factors, except H, become 
nonperturbative. 



B. Alternative Regulator, Consistent Subtraction and Relation to SCET Factor- 
ization 

In defining various contributions in Sterman's approach, the gauge choice A z = is made, 
or equivalently gauge links along the ^-direction are added to operators to make them gauge 
invariant. This choice of a non-light-like gauge can serve in addition gulator for 

collinear divergences arising from gauge links going along the light-cone direction, as can be 
seen from the one-loop the soft factor, Eq. (12*71) . 

In this subsection we present an alternative method to arrive at the correct factorization 
formula with factors that are manifestly gauge invariant. We regulate collinear divergences 



10 



by choosi ng gaug e links slightly off the light-cone (for more discussion of off-light-cone gauge 
links see |l3l E3. Il8|.) The direction of a gauge link supports a finite rapidity which can 
serve as a rapidity cutoff, thereby avoiding light-cone singularities which appear in calcula- 
tions of scaleless quantities like the parton distribution and which cannot be regularized by 
dimensional regularization We show that the factorization theorem proposed here is 
obtained only after proper subtractions of soft factors are made. By choosing the rapidity 
parameter p appropriately, we can eliminate the intermediate scale (1 — x)Q and arrive at 
a factorization similar to that of EFT. 

Let v = (v + ,v~,0) with v + ^> v~ and v = (v + ,v~,0) with v~ ^> v + with p = v + /v~ = 
v~ /v + . It is assumed below that the incoming quark is collinear in the z direction with 
momentum p\ = (Q/y/2, 0, 0) and the outgoing quark is collinear in the —z direction with 
p 2 = (0, Q/\/2, 0), and we denote p + = p{ = p 2 = Q/y/2. We define the soft factor as 

S{l-x) = [ % a ^ p+ (0\Ti[Y^0,-oo;Xv)Y^oo,0;Xv) 
J lix 

xY v (oo,0;0)y c (0,-oo;0)]|0)-^ (29) 

where Y v is 

Y v (a,b]£)=Vexp(-ig £ dX'v ■ A(X'v + m , (30) 

Similar definition holds for Y$. Thus the soft factor depends on two off-light-cone Wilson lines 
in the directions of v and v. This definition of the soft factor has no collinear divergences. On 
the other hand, if we take one of the Wilson lines on the light-cone, the resulting light-cone 
divergence may be considered as the collinear divergence. Then the S factor will include a 
collinear contribution as discussed in 3- However, here we are interested in the soft factor 



that is not contaminated with collinear divergences. 
The final state jet function is defined as follows: 

J(x,Q) = J ^e^-^^lW^Coo^jAuJV'CAu^W'^oo^jO)^) , (31) 

where it involves a Wilson line in the v direction which is taken to be in the (almost) 
conjugate direction to the out-going partons. It is given by: 

W v (a,b;C) =Pexp^-i# dX'v ■ A{X'v + f) J . (32) 

Finally we define the parton distribution function as 

m = A / ^e^ + (Pm^)W^(oo,0;Xv) 
2p + J 2tt 

x 7 + iy,(oo,0;0)^(0)|P) , (33) 

where the Wilson line is taken along the (almost) conjugate direction of the incoming partons, 
i.e., in the v direction. Both the jet and parton distribution are in principle defined to absorb 
just the collinear gluon contributions. Although this can be done by cutoffs in loop integrals, 
it is difficult to achieve in an operator approach. In fact, it will be clear later that both J(x) 
and 4>{x) do contain soft contributions as well, which must be subtracted explicitly. 
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Let us consider the factorization of the DIS structure function at one-loop using the 
above definitions of the factors. We first calculate the soft factor, jet function, and parton 
distribution in pQCD. The one-loop soft factor is 



S{l-x) = 5(1 - x) + ^C F ( -2 + ^llnp 2 ) (2D + hx%S{l-x 



2vr ~ * V ' p 2 - 1 



5(1 - x) + ^C F (-2 + hip 2 ) (2D + In ^5(1 - x^j 



(34) 



where the second line is obtained by taking the large p 2 limit. The result does not have any 
soft and collinear divergences. It does have an ultraviolet divergence coming from the cusp 
of the Wilson lines, which has been subtracted minimally in dimensional regularization. In 
moment space, 



^ = l + |^[2-lnp 2 ]ln^- 



(35) 



The UV-subtracted soft factor obeys the following renormalization group equation (RGE) 

dS(l -x,p 2 ^ 



dp 



2 ls S(l-x } p 2 ) 



where the anomalous dimension is 



ls = ^C F (2-\np 2 
In 



(36) 



(37) 



which depends on the rapidity cutoff p and is related to the so-called cusp anomalous di- 
mension 20] . 

The jet function has no infrared divergences either. At one-loop, 



a. 



J(x,Q) = 5(l-x) + -^C F 

Air 



3 - 2tt 2 - 3 In ^ + 4 lnp - 2 In 2 p ) 5(1 - x) (38) 



-7L> -4Di +41npA) 
In moment space, it has a particularly simple form, 



47T 



Np 



J N = 1 + -^C F { 31n-^ - 2[\n(Np) - l] 2 + 5 



Q 



-it 



(39) 



The wave function renormalization brings in the scale-dependence of the jet function, which 
therefore obeys the following RGE, 



/'■ 



dJ 



N 



dp 



with 



3a, 



7J,i 



2tt 

which is the anomalous dimension of the quark field in axial gauge 



(40) 
(41) 



21 
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The parton distribution at one-loop is, 



a* 



<t>(x) = 5(1 -x) + -^C F 
An 



Q 



4 + 



2tt 5 



6(1 -x) 



( 41n^ + 41np-4 ) D + 8A 



where the 1/e pole comes from collinear divergences. In moment space we have 



1 a * 
^ = 47 



3-41niV +4 



— — — — O 2 

■l+lniV + ln 2 ]V-lniVlnp-lniVln-\ 

A* 



(42) 



(43) 



Its UV divergences come from wave function corrections and have been subtracted minimally, 
and therefore its evolution in p 2 is the same as that for the jet function. 

Because and J contain the soft contribution as well, the structure function F\ cannot 
be factorized into <j> <g> J <g) S multiplied by a hard contribution, where <g> is a convolution 
operator in x-space. Instead, one must define the soft-subtracted version of the parton 
distribution <p = <p/S and jet J = J/S. Then the factorization reads: 



F ± (x,Q 2 ) = H(Q 2 )®J(x,Q 2 /p 2 



(x, Q 2 /p 2 ) ® S(x, Q 2 /p 2 ) + 0(1 - x) 



= H(Q 2 ) <g> J(x, Q 2 /p 2 ) ® 4>(x, Q 2 /fi 2 )/S(x, Q 2 /fi 2 ) + 0(1 - x) , (44) 

where H(Q 2 ) is hard contribution independent of x. This can be easily checked at one-loop 
level by calculating H, 



H W(Q 2 ) = Fi 1 \x,Q 2 )- J^(x,Q 2 )-^(x,Q 2 ) + S^(x,Q 2 ) 



1 + 



An 



o 2 o 2 

2 In \ + A In \ In p + 2 In 2 p - A In p - 8 + 2%' 

p l p l 



5(1 - x) (45) 



which is indeed independent of Di(x). All singular contributions of type 1/(1— x) + have been 
subtracted from the structure function F\ by the jet, soft factor and parton distribution. 
We emphasize that this is possible only when the soft contributions to the jet and parton 
distributions have been subtracted first. The anomalous dimension of the hard part is 



_ pdH a s , v 

1h = -T7-T- = C F (21np + 1) , 

H dp n ' 



(46) 



which is p-dependent. 

The above factorization uses the rapidity cutoff parameter p, which has the similar role 
as the renormalization scale p 2 : Every factor is a function of it, but the product has no 
dependence. Therefore, one can get different pictures of refactorization by choosing different 
value of p. For instance, if one takes p — > oo, all gauge links move back to the light- 
cone. Here collinear divergences shows up in different factors which have to be subtracted 
beforehand to yield a meaningful factorization. Distribution <fi corresponds to the physical 
quark distribution q(x). The subtraction of the soft contribution in J(x) ensures J(x) have 
collinear contributions only. Thus factorization can be written as 



Fi(x, Q 2 ) = H(Q 2 ) ® J(x, Q 2 /p 2 ) ® fa, Q 2 /p 2 ) + 0(1 - x) 
which is heuristically similar to Eq. (15) 



(47) 
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One can also take p = — 1. Then j(x) = J s t(x), and (j)(x) = 4> s t(x) and S(x) = (x), 
where quantities with subscripts "st" refer to those in the previous subsection. Eq. (44) 
then reproduces the factorization Eq. (20) from Ref. Q. 

We can also make contact with the EFT approach by taking p to be small although it 
shall be considered as large in principle. Consider the moment of <f)(x) = <p(x)/S, 



1 + ^C F {-1 [3 - AlnN] + 4 [l - In(iVp)] In - 4 } 



(48) 



where the finite part depends on a single logarithm In -^2, i.e., the scale (1 — x)Q. In the 

leading-logarithmic approximation, if p is taken to be 1/N, the distribution <j)(x) no longer 
depends on scale (1 — x)Q. The anomalous dimension of the distribution 0(x) is 

With p ~ 1/N, the evolution equation for <f)(x) is similar to that of the light-cone quark 
distribution q(x). 

Now let us examine the refactorization in the following form, 

F \{Q 2 ) = H(Q 2 /p 2 , p) ® J N (Q 2 /p 2 , p) ® <p N (Q 2 /p 2 , p) \ p „ 1/N + 0(1 - x) , (50) 

Taking p 2 = p 2 = Q 2 /N, the jet factor J N in Eq. (39) does not contain any large loga- 
rithms. The hard factor H(Q 2 / p 2 ) contains large logarithms that can be resummed. The 
resummation generates exactly the evolution of the matching coefficient C 2 for the product 
of effective currents in SCET. Therefore, the above form of factorization exactly reproduces 
the SCET result in Eq. (15). 



V. SUMMARY 

In this paper, we considered deep- inelastic scattering in a region (1 — x)Q ~ Aqcd where 
the standard pQCD factorization is not supposed to work. We argued, however, there is 
nothing that invalidates it in the new regime in leading order 1 — x because the Lorentz 
invariant factorization does not involve this soft scale in the sense that there are no new 
infrared divergences associating with this scale. 

We then discussed refactorization of the coefficient function. The EFT approach main- 
tains Lorentz invariance and hence allows a form of refactorization which is valid in the new 
regime. However, in the traditional approach in which jets and parton distributions are de- 
fined to take into account explicitly the double-logarithmic soft radiations, the scale (1 — x)Q 
does appear in various factors, making them nonperturbative in nature. We consider a more 
general factorization in this spirit which involves a rapidity cutoff. We showed how the EFT 
result can be reproduced through choices of this cutoff. The example also shows how to 
make consistent subtraction of the soft contribution in collinear matrix elements. 
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